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Recently spin-orbit (SO) coupled superfiuids in free space or harmonic traps have been exten- 
sively studied, motivated by the recent experimental realization of SO coupling for Bose-Einstein 
condensates (BEG). However, the rich physics of SO coupled BEG in optical lattices has been largely 
unexplored. In this paper, we show that in suitable parameter region the lowest Bloch state forms 
an isolated flat band in a one dimensional (ID) SO coupled optical lattice, which thus provides an 
experimentally feasible platform for exploring the recently celebrated topological flat band physics 
in lattice systems. We show that the flat band is preserved even with the mean field interaction in 
EEC. We investigate the superfluidity of the EEC in SO coupled lattices through dynamical and 
Landau stability analysis, and show that the BEG is stable on the whole flat band. 
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I. INTRODUCTION 



Flat bands possess macroscopic level degeneracy be- 
cause of their flat energy dispersion. They play a crucial 
role in important physical phenomena such as fractional 
quantum Hall effects where a large magnetic field applied 
on a two dimensional electron gas induces flat Landau 
levels [i|. The flat band physics is also greatly enriched 
recently by studying various lattice models where flat 
bands can be generated through geometrical frustration 
of hopping [2|, Q (e.g. Kagome lattice), the destructive 
interference between nearest-neighbor and higher-order 
tunnehngs (such as next-nearest-neighbor) [4-6], or the 
p-orbital physics [71]. In particular, isolated flat bands 
in lattices with non-trivial topological properties have 
attracted much attention in condensed matter physics 
for their applications in engineering fractional topologi- 
cal quantum insulator |8l-[l^ without Landau levels. 

However, most previous lattice models for generating 
fiat bands involve either high orbital bands or high order 
tunnehngs, which are generally very challenging in ex- 
periments. In this paper, we propose an experimentally 
feasible route for generating isolated flat bands using cold 
atoms in SO coupled weak optical lattices. Our work is 
motivated by the recent experimental realization of SO 
coupling for BEG [15], which opens a completely new av- 
enue for exploring SO coupled superfiuids [16]. In partic- 
ular, SO coupled BEG and degenerated fermi gases in free 
space and harmonic traps have been extensively investi- 
gated recently [l7l-l39|. However, ultra-cold atoms in SO 
coupled optical lattices have been largely unexplored [40] . 
We show that the combination of SO coupling, Zeeman 
field and optical lattice potential can yield isolated flat 
bands where topological properties may originate from 
the SO coupling [411, |42|. In regular optical lattices, the 



minimum of the lowest Bloch band locates at the center 
of the first Brillouin zone (BZ), while the maximum at 
the edge [43|. In SO coupled optical lattices, the mini- 
mum may locate at the edge and the peak at the center. 
The height of the central peak can be reduced with in- 
creasing Zeeman field, leading to decreasing band width 
and flat bands in certain parameter region. We note that 
such fiat band dispersion has been observed very recently 
in experiments in SO coupled optical lattices using ^Li 
Fermi atoms [44] . 

We first investigate a single atom in a ID SO coupled 
weak optical lattice to illustrate the mechanism for gen- 
erating isolated flat bands. The atom-atom interaction 
in BEG is then taken into account using the mean-field 
Gross-Pitaevskii (G-P) equation [45]. The nonlinear in- 
teraction reduces the band flatness, but does not fully 
destroy the isolated fiat bands. The combination of non- 
linear interaction and fiat bands may lead to rich and 
interesting physics. In particular, the instability of the 
nonlinear Bloch waves is very important because it di- 
rectly relates to the breakdown of superfluidity of the 
BEG [46l-[56j. In SO coupled optical lattices, the non- 
zero momentum of the energy minimum of the lowest 
Bloch band and the existence of fiat bands make their 
stability very different from regular optical lattices. For 
instance, the nonlinear Bloch waves can be stable in the 
whole BZ in the fiat band region. 

The rest of the manuscript is organized as follows. In 
Sec. II, we present the fiat band structure in SO cou- 
pled optical lattices. In Sec. HI, we discuss the effects 
of mean-field interactions and analyzes the stability of 
the BEG in SO coupled optical lattices. Sec. IV is the 
conclusion. 



II. FLAT BANDS IN SO COUPLED OPTICAL 
LATTICES 



* Corresponding Author, Email: chuanwei.zhang@utdallas.edu 



We consider a BEG confined in a ID optical lattice po- 
tential Vosin^ {kLx) along the x direction with Vq as the 




FIG. 1: (Color Online) Laser setup for implementing ID SO 
coupled optical lattices. Qi, and Q2 are the Rabi frequencies 
of the Raman lasers for generating SO coupling. The other 
two laser beams generate the optical lattice. 



lattice depth. In experiments, the lattice potential can 
be created by a standing wave formed by two lasers prop- 
agating along different directions [52] (see Fig. [1]). The 
effective wavevector of the lattice kL = 27rsin(6>L/2)/AL, 
where A^ is the wavelength of the lasers and Ol is the 
angle between two lasers. The SO coupling for BEC 
has been realized in experiments using two counter- 
propagating Raman lasers p^], yielding the single par- 
ticle Hamiltonian 






r^cr^. 



(1) 



where p is the atom momentum along the x direction, 
and a is the Pauli matrix. The SO coupling strength 
7 = hkR/m with kR = 27rsin(^i^/2)/Ai^, A^^ is the wave- 
length of the Raman lasers, and Or is the angle between 
Raman beams. Q is the Rabi frequency and acts as a 
Zeeman field. The units of the energy and length are 
chosen as the recoil energy 2El = h^k\/m^ and 1//cl 
respectively for the numerical calculation. Under these 
units, the single-particle Hamiltonian is dimensionless 
with 7 = kR/kL = sin(^i^/2)AL/sin(^L/2)Ai^ and the 
optical lattice potential l^ sin (x). 

Without optical lattice potentials, the single particle 
Hamiltonian Hq has two SO energy bands fi±{k) (shown 
in Fig. [2^) due to the lift of the spin degeneracy by 
the SO and Zeeman field. A gap 2Q between these two 
bands is opened at /c = by the Zeeman field Q. In 
the lower band, there are two energy minima at kmin = 
±^7^ — r^^/7^ and one peak at /c = 0. With increasing 
17, the distance between two kmin shrinks, and the height 
of the central peak decreases. At a critical value Qc = 7^ 
and beyond, two kmin merge to one point at kmin = 0, 
and the central peak vanishes. 

In the presence of periodic lattice potentials (i.e., con- 
sider the Hamiltonian Hq + Vo sin^(x)), the eigenenergies 
of the sing le-particle Hamiltonian form the Bloch energy 
bands [43|. To generate an isolated flat band, it is neces- 
sary to reduce the energy at both the edge and the center 
of the first BZ with respect to the band minimum, which 
can be realized through a combination of SO coupling, 
Zeeman field and lattice potential. Specifically, the pe- 




FIG. 2: (Color Online) Illustration of the formation of isolated 
flat bands, (a) Energy dispersion /j± (k) of the single-particle 
Hamiltonian Hq. The vertical solid (or dashed) lines corre- 
spond to the edges of the first BZ. (b) Lowest band for 7 < 1. 
The gray arrow indicates the suppression of the central peak 
with increasing Q. (c) The formation of the lowest band for 
7 > 1. The solid (green) line is /i_(/c), while the dotted lines 
are /i_(/c ± 2). 



riodic lattice potential can open an energy gap at the 
edge of the first BZ, which lowers the energy difference 
between the edge and the band minimum (denoted as h). 
When the original band minimum is close to the edge, the 
band edge becomes the band minimum (i.e., h = 0). On 
the other hand, the height of the central peak decreases 
with increasing Q. The band width should be determined 
by the larger value of h and the central peak height. 

The flat band generation mechanism is slightly differ- 
ent in two different regions: 7 < 1 and 7 > 1. For 
7 < 1, kmin = 7 at (7 = is within the first BZ (Fig. 
[JJd). If 7 is close to 1, /i should be zero and the band 
width is determined by the central peak height, which 
can be greatly reduced with increasing Q. Therefore the 
lowest band could be very flat for certain parameter re- 
gion. However, if 7 is much smaller than 1, h becomes 
a large value, and the width of the lowest band cannot 
be squeezed to the flat region. For 7 > 1, kmin of Hq 
lays outside of the first BZ. In this case, the lowest band 
is formed through folding the energy spectrum into the 
first BZ (i.e., shift the energy band of Hq by a lattice 
vector /c ± 2, see Fig. [2t). The band minima now locate 
at kmin T 2, and the physics is similar as that in 7 < 1. 
However, there is one major difference between 7 > 1 and 
7 < 1 . For 7 > 1 , the minimum of the lowest band first 
shift towards the edge of the first BZ when Q increases 
from 0. Therefore at certain range of Q, the band min- 
imum always stays at the band edge and the fiat band 
can be realized by suppressing the central peak with in- 
creasing ft. We emphasize that the resulting flat band is 
the ground state of the SO coupled lattice, which further 
enhances its experimental feasibility because atoms are 
usually adiabatically loaded to the lowest band in exper- 
iments [52]. Such SO mechanism for flat bands is very 
different from previous schemes in literature using high 
order tunneling or high orbital physics. 

The above intuitive physical picture agrees well 
with the numerical results. Using the Bloch theo- 
rem, the Bloch waves can be written as ^(x,t) = 
^{x) exp(— i/i (k) t-\-ikx)^ where ^{x) is the periodic part 
of the Bloch wavefunction, and /i (/c) is the eigenenergy. 
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FIG. 3: (Color Online) Flatness (ratio K) of the lowest Bloch 
band. Fo = 1- (a) 7 = 0.74 and c = 0. (b) 7 = 1.05. 
Solid line: c = 0; Dashed line: c = 0.05. The inset is a 
typical example of flat band in the nonlinear Bloch spectrum 
with Vl — 1.15, c — 0.05. (c) Plot of the dependence of the 
maximum flatness on 7 with c = 0. (d) The dependence of 
the maximum flatness on c with 7 = 1.05. 



which can be calculated using the standard central equa- 
tion. We measure the flatness of the lowest Bloch band 
by the ratio R of the gap between the lowest and first ex- 
cited bands to the width of the lowest band. In Fig. [3l we 
plot the flatness R with respect to ^ for two different 7. 
In the calculation, we use \r = 804.1 nm and Xl = 840 
nm which are typical for ^^Rb atoms in experiments |15| . 
The optical lattice potential is weak Vq = 2El to make 
sure the flat band does not come from the high lattice 
potential. For simplicity we choose 9l = tt, and consider 
two different Or: Or = tt corresponds to 7 = 0.74 (in Fig. 
[3^1), and Or = 7r/2 corresponds to 7 = 1.05 (in Fig. [SJd). 
We see the maximum flatness can reach nearly 20/1 for 
7 = 0.74 and 170/1 for 7 = 1.05. The suppression of the 
flatness for 7 = 0.74 < 1 agrees with our intuitive phys- 
ical picture: the band minimum for 7 = 0.74 is a little 
bit far from the edge of the first BZ, therefore the lowest 
band cannot be squeezed to exactly flat. In experiments, 
7 can be varied using laser setups with different Ol and 
Or or through a fast modulation of the laser intensities of 
the Raman lasers [32] . The dependence of the maximum 
flatness on the SO coupling 7 is plotted in Fig. [3b. With 
increasing SO coupling, the lowest band becomes more 
fiat. 



III. STABILITY OF EEC IN SO COUPLED 
OPTICAL LATTICES 

So far the study has been limited to the linear case, 
i.e., a single atom, while the interactions between atoms 
in BEC may play a major role on the dynamics of BEC. 



In the presence of a weak lattice potential, the mean 
field theory still applies and the dynamics of BEC in 
SO coupled optical lattices can be described by the non- 
hnear G-P equation 



5^ 



Ho^^Vosm\x)^^c{\^i\ 



|^2P)^, 



(2) 



where ^ = (^1,^2)"^ is the two component wave function 
of the BEC. The unit of time is m/hk\ and the wave- 
function is normalized through j dx{\^i\^ + l^2p) = 1 
in one unit cell. The dimensionless interaction coefficient 
c = h^uOyUJzkLdN / El^ where N is the atom number in 
one unit cell, a is the s-wave scattering length and uOy and 
ujz are the trapping frequencies in the transverse direc- 
tions. We consider a ID BEC confined in an elongated 
cigar-shaped trap with high transverse trapping frequen- 
cies {y and z directions), while the trapping potential in 
the longitudinal direction (x) is negligible. We also as- 
sume the interaction coefficients between atoms are the 
same for different hyperfine states, which is a very good 
approximation because their difference is very small [531 • 

Even in the presence of nonlinear terms, the solution 
of the GPE is still the Bloch wave in a periodic opti- 
cal lattice J46|, [43 • The repulsive interaction shifts the 
Bloch spectrum upwards, and modifies each band disper- 
sion and energy gap at the same time. However, isolated 
flat bands still exist in the presence of nonlinearity, as 
shown in Fig. [Sb where the flatness of the lowest band 
is plotted for c = 0.05. Compared with the linear case 
c = 0, the flatness of the nonlinear flat band decreases 
with increasing nonlinearity (Fig. [Sjl) and the maximum 
flatness is shifted towards a smaller value of 17. A typical 
example of the nonlinear Bloch spectrum with an isolated 
flat band is shown in the inset of Fig. [SJd. 

The combination of nonlinear interaction and flat 
bands may lead to various important phenomena, one of 
which is the superfluidity of the BEC in SO coupled op- 
tical lattices. For BEC in optical lattices, the breakdown 
of superfluidity may be caused by two different types of 
instabilities of the BEC, dynamical instability and Lan- 
dau instability, both of which have been extensively stud- 
ied in theory and experiments |46l-[56|. The existence 
of SO coupling and flat bands may significantly modify 
the superfluidity of the BEC. The stability analysis can 
be performed through Bogoliubov theory, where quasi- 
particle excitations induced by perturbations are taken 
into account through a small modification of the wave- 
function ^i(x,t) = [^i{x) -\-A^i{x,t)]ex.p{—ifit-\-ikx), 
where ^i{x) is the ground state of the BEC, A<l>^(x,t) = 
Ui{x) ex.-p{iqx — i5t) -\-w^{x) ex.-p{—iqx-\-iS^t)^ q and 5 are 
the wavevector and energy of the quasiparticle excita- 
tions, while Ui and Wi are the quasiparticle amplitudes. 
Substituting the modified wavefunction into the GPE, 
and linearizing the GPE with respect to Ui and Wi^ we 
obtain Bogoliubov-de Gennes (BdG) equations Sep = A4(p 
with (p = {ui,Wi,U2^W2)^ , where the matrix 



^~ \B2i A21 



(3) 
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FIG. 4: (Color Online) Instability of the lowest Bloch band, 
(a) and (b): Plots of the maximum of Im(^) for dynamical 
instability (a) and the negative maximum of /3 for Landau 
instability (b). Vo = 1, c = 0.01, 7 = 1.05, and Q = 0.8. 
(cl), (c2), (dl) and (d2): instability regions in the (/c, Q) 
plane. Dashed lines: the minimum of the lowest band. Blue 
grids: dynamical instability regions. Gray shadows: Landau 
instability regions. 
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Here C^'^^'^k) = -l/2{d/dx + ik)'^ + Vbsin^(x) - 
ij{d/dx-\-ik) — jj.-\-2c\^rn\'^ -\-c\^ri\'^. Because the matrix 
Ai is not Hermitian, its eigenvalues may be imaginary. 
The dynamical instability is defined if A4 has one or more 
non-zero imaginary eigenvalues. In this case, the insta- 
bility is characterized by the exponential growth of the 
perturbation. The nonlinear Bloch wave is dynamically 
stable if all eigenvalues are real numbers. On the other 
hand, the Landau instability can be studied by solving 
the BdG equation [46], Pep = TzM.^ with r^ = I^a^. The 
nonlinear Bloch wave ^i is said to be Landau instable if 
one or more eigenvalues of TzM are negative. Physically, 
the nonlinear Bloch wave with Landau instability is not 
the local energy minimum of the system. 

We systematically study the stability of nonlinear 
Bloch waves at the lowest Bloch band for various param- 
eters. A typical example of the dynamical and Landau 
instability of the lowest band is shown in Figs. H^ andHjD. 



Due to the symmetry in the plane (^, /c), we only show 
the region < (7 < 1 and < /c < 1. In Fig. |4^, the max- 
imum of the imaginary part of 3 is shown. The BEG in 
the region with non-zero values is dynamically unstable. 
There is a critical kd beyond which Bloch waves become 
dynamically stable. In Fig. HJd, the negative maximum of 
P is plotted, and non-zero values indicate the Landau in- 
stability. There is also a critical kc2 beyond which Bloch 
waves are the local energy minimum. 

In Figs. Ht and mi, we plot the dynamical and Landau 
instability region for different nonlinearity, SO coupling 
and Zeeman field. For BEG in a regular optical lattice, 
the energy minimum of the lowest band locates at /c = 0, 
and the Bloch waves in the region around /c = are 
stable. While in SO coupled optical lattices, the energy 
minimum of the lowest band may not locate at /c = 0, 
therefore we expect the stability domains should change 
accordingly, as clearly shown in Figs. Ht and |4]i. For 
7 = 0.74 < 1 in Figs. Htl and lit 2, the minimum of the 
lowest band (dashed lines) shrinks to A: = with increas- 
ing r^, and the abrupt change of the energy minimum 
corresponds to the flat band region. When 7 is close to 
one half of the first BZ, e.g., 7 = 1.05 in Figs. [Uil and 
[1112, the energy minimum initially stays at the edge of 
the first BZ k = 1. In the fiat band region, the energy 
minimum quickly moves to /c = 0. For a larger 7, the 
energy minimum initially increases from a value smaller 
than A: = 1 to the edge with increasing Vt^ stays there for 
certain range of Vt^ and then suddenly moves to /c = 0. 
The numerical results agree with the natural expectation 
that Bloch waves surrounding the minimum of the lowest 
band are stable, as shown in Fig. Ht and [Hi. However, 
we see the whole band is stable in the fiat band region, 
which means that the superfluidity of BEG with any mo- 
mentum in the fiat band is conserved. There are another 
two properties: 1) the region of dynamic instability is 
always smaller than the region of Landau instability; 2) 
the stable region increases for a larger nonlinear coeffi- 
cient c. These two properties are the same as those for 
BEG in regular lattices [ii,!!^. 



IV. CONCLUSION 

In summary, we show that the combination of SO cou- 
pling, Zeeman field and optical lattice can generate fiat 
ground state energy bands where the superfiuid of the 
BEG is stable in the whole band region. Our proposed 
SO coupling mechanism, when generalized to 2D, may 
provide an experimentally feasible route for generating 
chiral fiat bands and studying relevant fractional quan- 
tum Hall insulator physics. The stable superfluidity in 
the whole ground state band may lead to other interest- 
ing phenomena that have not been explored in regular 
optical lattices, such as dissipationless Bloch oscillation 
of BEG. 
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